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PREFACE 


This  Interim  Report  was  prepared  by  Wright  Aeronautical  Laboratories  at 
Wright-Patterson  Air  Force  Base  in  Ohio,  on  a  reimbursable  agreement  with  the 
U.S.  Department  of  Transportation,  Transportation  Systems  Center  (TSC) , 
Cambridge,  MA.  TAU  Corporation  of  Los  Gatos,  CA,  served  as  subcontractor  for 
this  project.  The  work  was  co-sponsored  by  the  U.S.  Coast  Guard  Office  of 
Navigation,  Office  of  Research  and  Development;  and  by  the  Research  and  Special 
Programs  Administration.  J.  H.  Kraemer  of  the  Center  for  Navigation  at  TSC  was 
the  Technical  Monitor,  assisted  by  R.  M.  Kalafus. 

The  algorithms  described  in  this  report  determine  the  differential 
corrections  to  be  broadcast  to  users  of  the (Global  Positioning  System  (GPSy  who 
require  higher  accuracy  navigation  or  position  information  than  the  30  to  100 
meters  that  GPS  normally  provides.  Properly  equipped  marine  users  receiving  the 
broadcasted  corrections  are  expected  to  achieve  better  than  10  meters  accuracy 
with  the  Differential  GPS  system. 

The  Differential  GPS  program  involves  the  fabrication  of  a  differential  GPS 
testbed,  consisting  of  an  8-channel  GPS  receiver  for  the  reference  station,  data 
processing  hardware  and  software  to  generate  the  corrections  and  perform  channel 
management  and  integrity  functions,  user  receiver,  and  data  acquistion  hardware 
and  software.  The  system  will  be  tested  in  1987. 


Other  reports  which  are  or  will  soon  be  available  include  a  specification 
of  the  design,  a  description  of  the  integrity  management  function  algorithms, 
and  test  results. 
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1.0  SYSTEM  CONCEPT 


The  differential  GPS  (DGPS)  reference  station  is  a  GPS  receiver  system 
fixed  to  a  known,  surveyed  location,  measuring  GPS  ranges  to  all  satellites  in 
view,  and  computing  and  broadcasting  range  bias  error  corrections  to  local 
users.  The  user  can  then  remove  the  locally  common  bias  errors  from  his  own 
measurements  before  computing  his  position/time  coordinates.  The  resulting 
DGPS  navigation  solution  is  substantially  more  accurate  than  a  simple  stand¬ 
alone  GPS  navigation  solution  with  only  minor  changes  to  the  user  equipment. 

2.0  ALGORITHM  REQUIREMENTS 

The  pseudorange  corrections  generated  by  the  Pseudorange  Correction  Gen¬ 
erator  (PRCG)  must  accurately  estimate  the  bias  errors  common  to  local 
users.  To  do  this,  the  reference  receiver  computes  nominal  range  to  the  sat¬ 
ellites  and  compares  against  the  measured  range.  Error  sources  that  are  not 
locally  common,  i.e.,  receiver-unique  error  sources,  must  be  removed  from  the 
measurements  before  broadcasting.  Table  1  contains  an  example  list  of  error 
sources,  both  receiver-unique  and  locally  common,  with  their  amplitudes,  fre¬ 
quencies,  and  correlation  distances.  Even  though  the  parameters  shown  in 
Table  1  are  only  rough  order  of  magnitude,  the  table  does  convey  the  relative 
comparisons  common  to  reference  receivers  and  supports  the  rationale  behind 
the  PRCG  algorithm. 


Table  1 .  GPS  Error  Sources 


Correlation 

Correlation 

Error  Source 

Amplitude 

Time 

Distance 

Receiver  Noise 

.01/10  m 

<10  sec 

0. 

( carrier/code) 
Interchannel  Bias 

1  ra 

500  sec 

0. 

Multipath 

20  m 

500  sec 

10  m 

Tropospheric  Delay 

30  m 

1000  sec 

20  km 

Satellite  Error  (SA) 

100  m 

100  sec 

400  km 

Ionospheric  Delay 

60  m 

1000  sec 

1000  km 

Satellite  Error  (no  SA) 

5  m 

10000  sec 

1000  km 

The  last  four  error  sources,  tropospheric  and  ionospheric  delay,  and  sat¬ 
ellite  error,  are  exactly  the  locally  common  bias  errors  which  DGPS  is  design¬ 
ed  to  eliminate.  (Satellite  error  includes  position/time  errors  resulting 
from  incorrect  satellite  orbit  and  clock  parameters  in  Space  Vehicle  (SV)  Data 
Blocks  I  and  II.)  The  other  three  error  sources  must  be  removed  before  the 
corrections  can  be  computed. 

The  first,  receiver  noise,  can  be  effectively  removed  by  filtering.  In 
fact,  the  heart  of  the  PRCG  algorithm  is  a  set  of  parallel  Kalman  (linear 
optimal)  filters,  one  filter  for  each  satellite  in  view. 

Interchannel  bias  is  f.  problem  unique  to  multichannel  receivers  such  as 
the  target  receiver,  the  Trimble  4000A.  The  magnitude  of  this  error  is  rela¬ 
tively  small  and  is  easily  measured  if  a  spare  channel  is  available.  The 
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spare  channel,  or  rover,  can  measure  itself  against  each  of  the  other  channels 
by  locking  in  on  the  same  satellite  and  comparing  measurements  over  a  period 
of  time.  The  average  difference  between  measurements  is  the  bias  relative  to 
the  roving  channel.  The  Trimble  4000A  performs  this  procedure  among  all  four 
channels  simultaneously  upon  initialization  and  can  also  be  commanded  to 
recalibrate  the  interchannel  biases  from  the  software  interface.  TAU  has 
studied  interchannel  bias  drift  in  the  4000A  and  has  determined  that  it  is 
negligible  over  a  six-hour  period.  For  the  demonstration  program,  no  inter¬ 
channel  bias  functions  will  be  implemented  other  than  recalibration  at  begin¬ 
ning  and  end  of  the  test  periods. 

The  -third  error  source,  multipath,  is  due  to  interference  from  reflected 
GPS  signals.  Multipath  effects  must  be  removed  since  they  ar“  generally  not 
common  to  the  reference  station  and  user.  However,  they  can  be  easily 
mistaken  for  ionospheric  delay,  so  they  must  be  removed  by  methods  other  than 
filtering.  The  most  effective  method  is  to  shield  the  antenna  from  all 
reflective  sources,  i.e.,  build  an  antenna  mask  that  allows  only  a  clear-sky 
view.  For  a  stationary  antenna  such  as  a  DGPS  reference  station,  this  is 
tractable.  Another  method  of  removing  multipath  effects  is  to  estimate  the 
multipath  error  and  compensate,  using  a  large  data  base  of  previous  days' 
readings.  This  may  be  possible  because  multipath  patterns  should  be  identical 
from  day  to  day,  assuming  the  satellite  rises/sets  in  the  same  direction  and 
the  multipath  sources  (buildings,  lakes,  etc.)  do  not  change.  The  third 
method  of  eliminating  multipath  is  a  multiple  antenna  approach,  assigning  a 
satellite  to  the  antenna  least  susceptible  to  multipath  for  that  satellite  at 
that  time.  The  current  demonstraton  program  assumes  multipath  has  been 
removed  by  the  first  method,  as  it  is  the  most  likely  to  be  implemented  in 
operational  systems. 

Upon  elimination  of  noise  and  multipath  effects,  only  locally  common 
errors  are  left.  Among  these,  troposphere  may  decorrelate  over  relatively 
short  distances,  so  special  provision  is  made  in  the  RTCM  format  for  distant 
users  to  recompute  a  tropospheric  correction.  Ionosphere  and  satellite  orbit 
errors  are  two  of  the  major  drivers  of  the  differential  concept,  because  they 
are  correlated  over  large  distances,  not  as  easily  modeled  as  troposphere,  and 
large  enough  to  preclude  precision  navigation. 

Table  1  reveals  a  major  requirement  for  proximity  in  a  temporal  sense. 
Most  error  sources  are  relatively  stable  over  time,  with  the  exception  of 
severe  Selective  Availability  conditions  and  possibly  rare  tropospheric  or 
ionospheric  disturbances.  These  cause  quick,  unpredictable  changes  In  the 
bias  errors  and  motivate  using  update  rates  that  are  higher  than  normal 
circumstances  would  dictate.  Without  these  unstable  error  sources,  error 
dynamics  are  fairly  small  over  a  few  minutes;  but  in  the  presence  of  these 
error  sources,  the  accuracy  of  the  corrections  may  degrade  over  seconds.  The 
solution  to  this  problem  requires  a  balance  between  expected  rate  of 
correction  degradation,  the  accuracy  requirements,  broadcast  message  content, 
and  available  bandwidth.  The  RTCM  format  compromises  on  a  nominal  50-bit  per 
second  data  rate  with  a  300-bit  message  (nominal  eight  satellites  in  view), 
yielding  a  six-second  update  rate.  Not  only  are  the  corrections  broadcast, 
but  also  the  rate  of  change  of  the  corrections,  so  that  temporal  degradation 
is  confined  to  the  second-order  effects  over  a  few  seconds.  Figure  1  shows 
correction  accuracy  degradation  for  this  scheme  over  time,  assuming  Selective 
Availability  behavior  matches  the  limited  data  sets  currently  available  to  the 
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public  [1],  A  six-second  update  rate  corresponds  to  accuracy  on  the  order  of 
five  meters,  sufficient  for  this  demonstration  project.  Six  seconds  also 
corresponds  neatly  to  current  GPS  broadcast  formats:  the  RTCM  message 
boundaries  coincide  with  the  Z-count.  This  provides  a  convenient  time 
reference  for  the  corrections  as  well;  the  broadcast  correction  is  the 
estimated  pseudorange  correction  for  the  particular  satellite  at  the  beginning 
boundary  of  the  message,  allowing  the  user  to  easily  extrapolate  it  to  the 
current  time  using  the  correction  rate  of  change  that  is  also  broadcast. 
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Figure  1.  Differential  GPS  Error  Growth  Due  to  Selective  Availability 

3.0  ALGORITHM  DESCRIPTION 

Each  satellite  in  view  has  a  dedicated  Kalman  f liter/ correction  genera¬ 
tor.  As  illustrated  in  Figure  2,  the  only  co-mingling  of  satellite  data 
occurs  in  the  clock  correction  filter,  which  is  a  heavily-damped  average  of 
all  pseudorange  measurements.  The  clock  correction  filter  is  described  in 
more  detail  after  the  following  sections  describing  the  individual  Kalman 
filters. 

The  Pseudorange  Correction  Generator  accepts  as  inputs  pseudorange  and 
deltarange  measurements,  satellite  ephemeris  data  and  clock  correction 
estimates;  and  outputs  pseudorange  error  estimates  (to  the  clock  filter), 
pseudorange  correction  and  correction  rates  of  change,  and  difference  between 
old  and  new  ephemeris.  Residuals  from  the  filter  are  stored  also  and  used  for 
integrity  management  purposes. 
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Figure  2.  Pseudorange  Correction  Generation 
3.1  Measurement  Processing 

The  raw  measurements  from  the  receiver  are  pseudorange  and  deltarange; 
before  input  to  the  filter,  they  must  be  transformed  to  pseudorange  and  delta- 
range  error  by  subtracting  the  nominal  range  and  change  in  range  as  well  as 
the  SV  clock  corrections.  The  data  to  perform  these  computations  come  from  SV 
Data  Blocks  I  and  II;  the  general  data  flow  is  illustrated  in  Figure  3.  The 
formulas  for  clock  correction,  satellite  position,  and  range  are  in  Tables  2, 
3,  and  4,  respectively. 

Nominal  deltaranges  are  obtained  by  computing  the  nominal  pseudoranges  at 
the  beginning  and  the  end  of  the  averaging  time  interval,  and  then  backwards 
differencing. 


Figure  3.  Pseudorange  Computation 


Table  2.  Satellite  Clock  Corrector 


Variables 

tot ,  tor 


SV,  RVR,  GPS 
P\ieas 


At 


sv 


toe,  afQ,  afj,  af2, 
e,  sin  E 
Equations 


tot 


SV 


At 


SV 


t°CGPS 


-  time  of  transmission,  reception 

-  SV,  receiver,  GPS  time  frames  (subscripts) 

-  measured  pseudorange 

-  SV  clock  error 

tgd  -  SV  Data  Block  I  parameters 

-  orbital  parameters  (cf.  Table  3) 

-  totRVR  -  PRmeas 

■  (afQ+tgd)  +  af  jCtotgy-toc)  +  af  2(totgV-toc)2 

+  (-4.42809305  x  10-10)  e  /a  sin  E 

-  totsv  -  Atgv 


I 
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Table  3.  Satellite  Orbit  Equations 


Variables 


u.n 


toe,  M  ,  A,  An,  e,  <»,  C  , 


^uc1  ^rs*  ^rc'  ^is»  ^ic’ 

i  ii  ij  a 

o».dt’  o’ 

M,  E,  v 


u,  r,  1,  ft 

xsv»  Ysv»  zsv 


-  earth  gravity,  rotation  constants  (WGS  72) 

-  SV  Data  Block  II  parameters 


-  mean,  eccentric,  true  anomalies 

-  inertial  orbital  elements 

SV  position  in  ECEF  at  time  of 
transmission 


Equations 

t 

M 

M 


totGps  -  toe 

M0  +  +  An)  t 

A 

E  -  e  sin  E 


tan  1 
v  +  w 


/(~I-e2)  sin  E 
cos  E  -  e 


solved  by  Newton's  root-finding  method, 
precision  <10-12  radians 

} 


u  =  <fr  +  C  sin  2d>  +  C  cos  2$ 
us  uc 

r  »  A(I-e  cos  E)  +  C  sin  2d>  +  C  cos  2$ 

rs  rc 

i  =  i„  +  (44)  t  +  C,  sin  2 $  +  C.  cos  2$ 
o  dt  is  ic 

a  =  a  +  (S  -  n  )t  -  5  t 

o'  e  e  ce 

Xgy  =  r  cos  u  cos  ft  -  r  sin  u  cos  i  sin  ft 

ygV  »  r  cos  u  sin  ft  +  r  sin  u  cos  i  cos  ft 

2gV  =  r  sin  u  sin  i 


6 


1 


Table  4.  Nominal  Pseudorange  Computation 


Variables 


SV,  RR 


XTX ’  XRX 

e 

r 

prnom 

rto 

Equations 


earth  rotation  rate,  speed  of  light 
ECEF  satellite,  receiver  positions 

vector  coordinates  in  ECEF  at  time  of  transmission  and 
reception 

angle  earth  rotates  between  t^,  t^ 
range  between  3V  and  reference  receiver 
nominal  pseudorange 

receiver  time  offset  (dual  receivers  only) 

SRrxI  •  Vc 


SV 


RX 


cose  sine 
-sine  cose 
0  0 


SV, 


TX 


SVRX  -  "‘Sot1 


prnom  *  r  +  Atsv  "  rt° 
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Note  that  nominal  pseudorange  as  defined  here  does  not  include  the 
effects  of  ionosphere  or  troposphere,  nor  does  it  include  the  effect  of  recei¬ 
ver  clock  error.  Therefore,  these  effects  will  be  included  in  the  pseudo- 
range/deltarange  error  measurements  entering  the  filter.  The  clock  error 
alone  is  on  the  order  of  hundreds  of  kilometers,  dwarfing  all  other  pseudo¬ 
range  error  sources  combined.  As  a  result,  the  clock  error  is  common  to  all 
satellites  in  view,  however,  and  is  estimated  and  removed  after  filtering. 


3 . 2  Kalman  Filter  Processing 

Each  satellite  in  view  is  assigned  a  dedicated  four-state  Kalman  fil¬ 
ter.  ThA  filter  processes  pseudorange  and  deltarange  measurement  upon  receipt 
and  outputs  a  pseudorange  error  function,  quadratic  in  time,  suitable  for 
input  to  the  pseudorange  correction  formulas. 


The  filter  will  be  defined  here  according  to  the  Kalman  filter  equations 
shown  in  Table  5  [2].  The  actual  arithmetic  implementation  is  the  U-D  formu¬ 
lation  laid  into  a  standard  software  version  by  G.  J.  Merman  [3].  Using  the 
terminology  of  Table  5,  the  filter  is  defined  as  follows: 


x^  ■  state  vector  at  time  tk  = 


PRE 

PRE 

PRE 

i 


where  PRE  -  pseudorange  error 

P&E  -  pseudorange  error  rate 


PRE  -  pseudorange  error  acceleration 
I  -  ionospheric  delay  rate 


w^  =  state  noise  tk_j  to  tk  = 


k-1 


4.  ,  =  transition  matrix 
k-1 

tk_!  to  tk  = 


0 

0 

PR  acceleration  noise 
ionospheric  delay  rate  noise 


1 

0 

0 

0 


^k^k-P 

1 

0 

0 


1 / 2 ( tk_tk-l )2  0 

(tk“ck-l)  0 
1  0 

0  1 


Qk  =»  state  noise  covariance  = 


k-1 


0 

0 

0 

0 


0  0  0 

0  0  0 

0  qp  0 

0  0  qx. 


with  qp 


(TBD)  and  qt  -  (TBD) ; 


^  =  (pseudorange 

error 

measurement,  deltarange  error  measurement); 

Hk  -  [l 

0 

0  o]  ; 

L° 

1 

0  -2J 

^k  "  pcode 

r  1 

L° 

rcarrierl 

with 

rcode 

(TBD)  and  rcarrier  *  (TBD). 
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Table  5.  Summary  of  Kalman  Filter  Equations 


System  Model 

Measurement  Model 

State  Estimate  Extrapolation 

Error  Covariance  Extrapolation 

State  Estimate  Update 

Error  Covariance  Update 

Kalman  Gain  Matrix 


5k  "  V  1-Vi  +  Sk-1  *  VN(2’V 
5k  =  \5k  +  VN<°-V 
5k<->  "  Vl  -V>(+) 

pk(->  -  ViVi«ViT  +  Vi 
ik(+)  =  5k(_)  +  V5k  " 
pk(+)  -  t1  -  \W’> 

Kk  =  pk(->  HkTt\pk(-)HkT  +  V'1 


The  Kalman  filter  has  three  unique  features  worthy  of  comment.  First, 
even  though  linear  corrections  are  being  broadcast,  quadratic  corrections  are 
being  computed.  As  will  be  demonstrated  in  the  section  on  pseudorange  correc¬ 
tion  computation,  this  significantly  lowers  the  latency  error.  This  latency 
error,  due  in  part  to  delays  in  transmission  and  in  part  to  the  six-second 
applicability  error,  becomes  less  sensitive  to  pseudorange  acceleration 
because  these  effects  are  estimated  and  partly  compensated. 

The  second  feature  is  the  presence  of  the  ionospheric  delay  rate  term. 
This  term  is  included  to  model  the  deltarange  measurement  properly.  As  shown 
in  Reference  4,  the  delta-range  is  not  truly  the  change  in  pseudorange: 

&PR  -  Ai  *  dr  +  Ai 

or  Apr  -  2A1  =  DR 

where  Apr,  Aj  are  change  in  pseudorange  and  ionospheric  delay,  respectively, 
and  DR  is  the  so-called  deltarange  as  measured  by  the  carrier  loops.  Although 
the  error  introduced  by  ignoring  ionospheric  delay  is  small,  it  can  build  up 
over  time  to  several  meters.  Rather  than  mismodel  the  system  with  no 
ionosphere  and  compensate  with  additional  process  noise,  the  extra  state  has 
been  added. 

The  last  feature  deals  with  Kalman  filter  process  noise  modeling  and  the 
character  of  the  measurements.  Typical  Kalman  filter  states  are  current 
Instantaneous  position,  velocity,  and  acceleration;  this  filter  has 
instantaneous  position  and  average  velocity  and  acceleration  over  the  time 
interval.  The  particular  formulation  of  the  state  is  dictated  by  the 
deltarange  measurement.  In  the  target  receiver,  the  deltarange  measurement  is 
derived  from  the  average  Doppler  over  the  interval  since  the  last  measurement, 
as  opposed  to  Instantaneous  Doppler  at  the  time  of  measurement.  Thus,  the 
deltaranges  reflect  the  total  change  in  range  over  time  instead  of  regularly 
sampled  range  rates.  For  position  estimation,  the  "average"  measurement  gives 
slightly  better  results  than  the  "instantaneous"  measurement  because  of 
slightly  lower  process  noise.  Consequently,  the  Trimble  4000A  performs  range 
tracking  slightly  better  than  receivers  using  an  instantaneous  range-rate  or 
Doppler. 
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3.3  Clock  Correction  Filter 


To  reduce  the  dynamic  range  of  the  pseudorange  corrections,  the  average 
of  the  pseudorange  errors  is  computed  and  removed  before  computing  the 
pseudorange  corrections.  Each  pseudorange  error  is  typically  the  sum  of  a 
very  large  common  term  (the  receiver  clock  error)  and  the  satellite-unique 
term  (comprising  ionospheric  and  tropospheric  delays  and  satellite  clock  and 
orbit  error).  Because  of  atmospheric  delays,  the  pseudorange  error  average 
tends  to  be  a  bit  larger  than  the  receiver  clock  error.  The  moniker  "clock 
correction"  is  used  because  the  difference  between  the  two  is  so  small  that 
correction  is  almost  exactly  equal  to  the  receiver  clock  error.  Any  reference 
receiver  clock  error  remaining  in  the  pseudorange  correction  translates  into 
an  equal  and  opposite  user  receiver  clock  error  when  the  differential 
corrections  are  applied.  Indeed,  simple  theoretical  principles  dictate  that 
the  correction  term  can  be  any  number,  as  long  as  it  is  subtracted  equally 
from  each  of  the  differential  corrections;  and  that  the  only  term  affected  is 
user  time,  not  user  position.  This  is  due  to  the  nature  of  differential  GPS, 
that  the  user  is  navigating  relative  to  the  receiver  in  position  as  well  as 
time,  so  a  receiver  clock  error  (or  error  in  the  receiver  clock  correction) 
results  in  only  a  user  time  error. 

In  practice,  the  above  is  nearly  true;  in  fact,  it  is  reliable  for  user 
receivers  such  as  the  Trimble  4000A,  which  make  no  effort  to  model  the  dyna¬ 
mics  of  the  user  clock.  For  more  sophisticated  receivers,  or  receivers  aided 
by  and  depending  on  cesium  clocks,  large  clock  error  dynamics  introduced  by 
the  differential  corrections  and  left  unmodeled  by  the  user  navigation  filter 
may  correlate  into  the  position  estimate  and  contaminate  the  navigation 
accuracy.  The  amount  of  contamination  depends  upon  how  heavily  the  user 
filter  relies  on  the  low  user  clock  error  dynamics. 


At  issue  in  the  design  of  the  reference  station  is  whether  the  reference 
receiver  clock  error  remaining  in  the  differential  corrections  has  high  enough 
unmodeled  dynamics  to  affect  the  user  position  solution.  The  most  reliable 
approach  is  to  keep  the  actual  clock  dynamics  low  with  a  rubidium  frequency 
standard,  as  is  being  done  in  this  project.  Once  the  clock  dynamics  have  been 
guaranteed  to  be  very  low,  a  pseudorange  average  will  be  able  to  follow  the 
true  reference  receiver  clock  error  fairly  closely  yet  not  contaminate  user 
position  error.  The  average  must  be  sufficiently  damped  so  that  impulses 
caused  by  rising  or  setting  satellites  or  uploads  cause  minimal  disturbance. 
A  simple,  two  state  filter  is  sufficient  to  create  this  damped  average.  The 
inputs  to  the  filter  are  the  simple  averages  of  the  current  pseudorange  and 
deltarange  errors.  A  suboptimal  constant  gain  linear  filter  will  be  used: 


a  a  «  - 

cc .  .  =  cc.  +  cc  .At  +  k.  PRE, 
1+1  11  11  l 


A 

cc . 


=  cc . +  k,_  PRE . 
i+l  i  22  i 


where 


ccj_  =  clock  correction  estimate  at  t^ 


At  ^  =  time  since  last  measurement  cycle  *  ti+l-ti 

PRE .  =  average  of  pseudorange  error  estimates  at  t^ 
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kll  =  k22  =  ~n“’  where  n  =  (TBD)* 

A  "corrected"  pseudorange  error  and  error  rate  are  then  generated: 


A 

A 

A 

PREc  - 

PRE 

-  CCi 

A 

A 

/Cv 

PREci  = 

PRE1 

-  CCj 

PREc .  =  PRE. 
l  i 

Acceleration  is  unaffected  because  rubidium  clock  frequency  drift  is  so  small. 


3.4  Computation  of  Linear  Pseudorange  Correction 

Each  quadratic  pseudorange  error  function  must  be  converted  to  a  linear 
pseudorange  correction  function.  Typically,  there  is  a  lag  between  the  Last 
available  data  input  into  the  filter  (denoted  t^)  and  the  start  of  the  broad¬ 
cast  interval  (denoted  t^)  in  which  the  correction  will  be  valid,  and  there 
sometimes  is  variation  in  the  length  of  the  time  interval  (At).  (See  Figure 
4.)  These  are  all  considered  in  the  following  algorithm,  which  computes  the 
line  of  best  fit  (in  a  maximum  error  sense)  to  the  quadratic  function  pre¬ 
dicted  by  the  pseudorange  error  filter. 


The  predicted  quadratic  coefficients  are  output  from  the  filter: 
PREc(t)  -  PREc1  +  PREc1(t-t1)  +  PMIc1(t-ti)2/2. 


The  rate  of  change  of  pseudorange  correction  is  equal  to  the  slope  of  the  tan¬ 
gent  of  the  PRE  curve  at  the  midpoint  of  the  broadcast  interval: 

.  A 

PRC  -  PREci  +  PREc1(tb-t1+At/2). 


Lastly,  the  value  of  PRC  is  the  ordinate  of  the  best  fit  line  at  time  tb.  The 
best  fit  line  is  midway  between  the  above-mentioned  tangent  and  the  secant 
joining  the  points  of  the  PRE  curve  at  times  t^ ,  tb+At: 

A  2 

PRC  -  PREc  (tb)  +  PREc1(At)  /16. 


f 


Figure  4.  Linear  Pseudorange  Correction 


11 


4.0  SUMMARY 


The  above  algorithm  takes  the  satellite  measurements  and  processes  them 
through  a  filter  to  remove  the  receiver-unique  error  sources,  which  tend  to  be 
higher  frequency  than  locally  common  error  sources.  An  estimate  of  clock 
error  is  computed  and  then  removed  equally  from  all  filter  outputs.  The 
clock-corrected  filter  outputs  are  coefficients  of  a  quadratic  function  in 
time  estimating  pseudorange  error.  A  best-fit  linear  approximation  to  this 
function  over  the  applicable  interval  is  generated;  the  coefficients  of  this 
linear  time  function  are  the  elements  to  be  broadcast  as  a  pseudorange 
correction. 
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